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ABSTRACT 

Using the publicly available VESPA database of SDSS Data Release 7 spectra, we 
calculate the stellar Mass Weighted Age (hereafter MWA) as a function of local galaxy 
density and dark matter halo mass. We compare our results with semi-analytic models 
from the public Millennium Simulation. We find that the stellar MWA has a large 
scatter which is inherent in the data and consistent with that seen in semi-analytic 
models. The stellar MWA is consistent with being independent (to first order) with 
local galaxy density, which is also seen in semi-analytic models. 

As a function of increasing dark matter halo mass (using the SDSS New York 
Value Added Group catalogues), we find that the average stellar MWA for member 
galaxies increases, which is again found in semi-analytic models. Furthermore we use 
public dark matter Mass Accretion History (MAH) code calibrated on simulations, 
to calculate the dark matter Mass Weighted Age as a function of dark matter halo 
mass. In agreement with earlier analyses, we find that the stellar MWA and the dark 
matter MWA are anti correlated for large mass halos, i.e, dark matter accretion does 
not seem to be the primary factor in determining when stellar mass was compiled. 
This effect can be described by down-sizing. 
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1 INTRODUCTION 

The spectra of galaxies encodes information about the his- 
tories of the component stellar populations, dust, and star 
formation. Various tools have been developed to extract this 
information (e.g.. Heavens et al. 2000 Tojeiro et al. 20071 



and previous works have compared the resulting extracted 
information with both extrinsic and intrinsic galaxy prop- 
erties. These approaches rely on the assumption that the 
evolution of the stellar populations are well understood and 
that the current modeling of stellar population is accurate. 

The MOPED ( [Heavens et al.|2000[ ) routine implements 
the general process of reforming a complex dataset (e.g., a 
galaxy spectra) into a set of parameters (e.g., star forma- 
tion rate, metallicity) and parameter combinations, assum- 
ing uncorrelated noise, such that the data compression is loss 



less (see also [Graff et al.||2010| ). [Mateus et al.| ( |2008| ) used 
MOPED-derived stellar masses and luminosities of SDSS 
Data Release 3 ( [Abazajian et aL][2005[ ) galaxies to build 
marked correlations functions, and compared with marked 
correlations of semi analytic models of galaxies in the Mil- 
lennium Simulation. 



Component Analysis technique, to decompose the spectra 
of low redshift {z < 0.1) SDSS early-type galaxies into two 
quantities: the average stellar age for the galaxy, assuming 
metal-rich stellar populations, and the fraction of recent star 
formation. They found little dependence (and a large scat- 
ter) of recent star formation and average stellar age on host 
halo mass, but find a stronger correlation on local processes, 
e.g., galactic velocity dispersion. 

An easily accessible, robust code, i s the VErsatile SPec- 
tral Analysi^ (hereafter VESPA, see Tojeiro et al 



2007 



[2009[ for more details) package, which recovers star forma- 
tion and metallicity histories of the galactic spectra using 
synthetic stellar population models. The software recovers 
histories in adaptive age bins according to the signal-to- 
noise of the galaxy spectrum on a case by case basis and ad- 
dresses the age-metallicity relation. Two popular synthetic 
stellar population models are included in the VESPA out- 
hereafter BC03), 



put, those of Bruzual & Chariot (2003 



and the Maraston ( 2005 I & Maraston et al. ( 2009 1 hereafter 



M05, which differ in their respective resolutions, and the 



More recently. 



Ferreras et al. 



(20101 applied a Principal 



^ http: / / www- wfau. roe .ac.uk/ vespa / 
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use of empirical libraries to model the thermally pulsating 
Eisymptotic giant branch. Furthermore VESPA corrects for 
galactic extinction using the dust maps of 'Schlcg el et al.| 
(19981, and fits for the dust in each galaxy using either a 



Tojeiro et al. ( 2009 1 to 



one or two parameter dust model. 

VESPA has been used by, e.g. 
compare the BC03 and M05 models, in particular, the re- 
covered total stellar mass today (finds good agreement) and 



the mass averaged metallicity (finding less agreement). To- 



jeiro et al. (20101 used VESPA to model the color evolution 



of high signal-to-noise SDSS Data Release 7 Luminous Red 
Galaxies as a function of color, luminosity and redshift. 

In this paper we compare the results of VESPA, in par- 
ticular the typical time scale of stellar mass assembly, the 
stellar Mass Weighted Age (hereafter MWA), with similar 
results from N-body simulations and semi analytic models. 
In particular we address the following questions: 

(i) Does the stellar MWA depend on the local density? 

(ii) Are similar trends seen in semi analytic models? 

To this end, we examine the correlation of local galaxy den- 
sity with the VESPA calculated stellar MWA, and compare 
with the stellar MWA of the semi analytic models in the 
Millennium Simulation ( [Springel et al.|2005||Boylan-Kolchin| 
et al.|2009| |. 



(iii) Is the stellar MWA correlated with the mass of the 
dark matter halo that the galaxy inhabits? 

In this case, we match the galaxies used in VESPA with 
those from the New York Value Added Group catalogue 



(Blanton & others. 2005 Yang et al. 2007) and explore the 



average stellar MWA of member galaxies as a function of 
dark matter halo mass. 

(iv) Is the stellar MWA correlated with the dark matter 
MWA? 

We address this by comparing the results of the above anal 
ysis with the results of Mass Accretion History code ( Zhao 



et al. 2009 1, rebinned in terms of VESPA time bins for direct 



comparison to the observations. 

The layout of the paper is as follows: in [|2]we describe 
the observational and simulated data and the process of 
reconstructing the stellar MWA of a galaxy. We continue 
in ^ by detailing our measurement of local galaxy den- 
sity and dark matter halo mass, for both the observed and 
semi analytic samples, and present the results in ij4] We 
conclude and discuss in i|5] Throughout the paper we em- 
ploy a flat ACDM cosmology with {h, JIa, fim, os) given by 
(0.7kms-iMpc, 0.7,0.3, 0.8) 



2 DATA 

In this paper we compare observational data with semi an- 
alytic models from simulations, and dark matter mass ac- 
cretion code calibrated on dark matter simulations, and we 
describe all data sets below. We also include the Structure 
Query Language (hereafter SQL) code to retrieve the data 
sets, and present the charateristics of the final data sets in 
Tabled 



2.1 Observed Galaxy Data 

The 6 X 10'' spectroscopically selected galaxies in this work 
were drawn from the Sloan Digital Sky Survey Data Re- 
lease 7 (York et al. 2000 Abazajian et al. 2009 and ref- 



erences therein, hereafter SDSS DR7). Galaxy targets were 
selected from SDSS imaging and include "Main Galaxy Sam- 
ple" (Strauss et al. 2002 hereafter MGS), and "Luminous 



Red Galaxies" ( Eisenstein et al.|2001 hereafter LRG) 



Galajcy spectra were reprocessed using VESPA and we 



followed Tojeiro et al. ( 2009 1 in adopting the two parameter 



dust model. We note that our results are insensitive to the 
choice of dust models but are more sensitive (to within 20%) 
to the choice of stellar population models, which we discuss 
in later sections. 

We obtained the galaxy data using the following query 
run on the VESPA SQL web interfac^ modifying the values 
of runid to obtain different dust models, stellar population 
models, and galaxy populations. We show the size of the 
galaxy samples obtained in Table [l] and describe the main 
parts of the query later in this section. 

SELECT tables. indexP, tables. AgeMass/table3.TotMass AS MassWeightedAge , 

tables. AgeZ/tableS.TotZ AS ZWeightedAge , 

tables. AgeSFR/tableS.TotSFR AS SFRWeightedAge , 

tables. TotSFR, tableS.TotZ, tables . TotMass , 

1 . specDbj ID, n, redshift ,d . dustval ,d. runid, d, dust id, 

n.m_stellar AS TotalCurrentStellarMass ,n .nbins , 

s.dered.g AS g, s .modelMagerr_g, s . dered_r AS r, s ,modelMagerr_r , 

s.dered_i AS i, s . modelMagerr_i , s.dered_z AS z, s.z AS spec_z, 

s . zErr AS errspec_z , s . modelMagerr_z , 

s.ra AS gal_ra, s.dec AS gal_dec, 

S.HJD, S.FIBERID, S.PLATE, S.TILE, z.absMag.r 

FROM ( 

SELECT indexP, sumCmass+meanAge) AS AgeMass, suinCmass) AS TotMass, 
sumCZ) AS TotZ, sumCZ*meanAge) AS AgeZ, suin(SFR) AS TotSFR, 
sumCSFR*raeanAge) AS AgeSFR 

FROM 

( SELECT b.indexp, b.mass AS mass, 

CbID . ageStart+bID . ageEnd) /2 AS meanAge, b.Z, b.SFR 

FROM binprop AS b 

JOIN binlD AS blD ON 

( b.binID - blD.binlD AND b. runid = 8) 

JOIN lookupTable AS 1 ON ( b.indexP - l.indexP) 

) AS table2 group by indexp 

) AS tables 

JOIN lookupTable AS 1 ON ( tables . indexP = l.indexP ) 

JOIN galprop n ON (tables . indexP=n. indexP AND n.runid=8) 

JOIN dustprop d ON CtableS . indexP=d. indexP AND d.runid=8 AND dustid=2) 

JOIN bestDR7. .SpecPhotoAll AS s ON Cs . specDbj ID=1 . specObj ID ) 

JOIN bestDR7. .photoz AS z ON (z . ObjID-s . ObjID ) 

In this work, we are concerned with one possible output 
of the VESPA database, namely the stellar Mass Weighted 
Age (MWA) . It is calculated, as in the above SQL query, as 
the sum over VESPA age bins i, of the mass Mi, formed in 
the bin, multiplied by the central bin age ti, divided by the 
total stellar mass formed in the galaxy i.e. 



MWA 



Mto, 



(1) 



The adaptive binning inherent to VESPA adopts wide 
age bins for spectra with low signal-to-noise, leading to an 

^ http://www-wfau.roe.ac.uk/vespa/SQL_form.htmI 
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VESPA runs on SDSS D7 galaxy spectra 



Description 


SP Model & runID 


Redshift range 


Number of galaxies 




BC03 2 


0.05 ^ 2 ^ 0.20 




T Rn C;amT-.lr. (1 P \ 
l^lVvj OclIlipiC (^ijIXVj J 


BC03 6 


0.35 ^ 2 ^ 0.53 




Main Galaxy Sample (MGS) 


M05 4 


0.05 ^ 2 ^ 0.20 


622749 


LRG Sample (LRG) 


M05 8 


0.35 ^ 2 ^ 0.53 


75423 


Semi Analytic 


galaxies from the Millennium Simulation DeLucia2006a 


table. 


Redshift 


Number of galaxies 


Number of parent DM Halos 




2 = 0.06 


21155 


7344 




z = 0.12 


22175 


8272 




z = 0.17 


21818 


8230 




z = 0.36 


138 


62 




2 = 0.41 


403 


138 




2 = 0.51 


479 


154 





Table 1. A description of the data samples used in this work. We show the acroymn used throughout this paper, the "runID" from 
the VESPA database applied to the SDSS DR7 spectra, the redshift cuts, and number of galaxies which have greater than 4 recovered 
VESPA age bins. We also show the galaxy data drawn from semi analytic models retrieved from the Millennium Simulation, the redshift 
of output simulation, the number of galaxies and the number of parent halos. 



unnatural peak in the MWA distribution at 7 Gyrs (which 
corresponds to the central value of the largest time bin, see 
Fig. 1 of Tojeiro et al. 20091. To remove these unphysical 



solutions we further imposed that the number of recovered 
VESPA bins is greater than 4, when determining averaged 
stellar MWA. We use the full data sets when calculating the 
local galaxy density. In Table [T] we present the final redshift 
ranges and the number of galaxies for both the MGS and the 
LRG samples using the different stellar population models. 



2.2 Simulated Galaxy Data 

The simulated data were drawn from the Millennium Sim- 
ulation, which used a modified version of the colhsionless 
dark matter and gas dynamical cosmological sim ulation code 
GADGET2 ( [Springel et al.|200l| |Springel|2005| with 2160^ 



dark matter particles, each with mass 8.6 X 10" /i"^Mq, m a 
volume of side 500 /i~^Mpc. The simulation is able to recover 
dark matter halos small enough to host luminous galaxies 
brighter than O.IL*. 

The dark matter halos were populated with galaxies ac- 



cording to the semi analytic models of De Lucia et al. ( 2006 1 



and De Lucia & Blaizot ( 2007 1. The models are novel in that 



they continue to follow dark matter halos until they are 
tidally disrupted by parent halos, and allow central galaxies 
to replenish their supply of cold gas, unlike satellite galaxies. 
They contain prescriptions for gas infall onto halos, mixing 
between hot and cold gaseous states and cooling to form 
stars, and allow for cooling flows and their suppression by 
central galaxy Active Galactic Nuclei activity (see Springel 
let al.||2005l [Croton et ar][2006l |De Lucia et al.||2004[ for 
specific details). 

Both the gala:xy properties and host dark-matter halo 
properties are stored in the "DeLucia2006a" table found on 
the public data access websit^ We obtained information 
about local galaxy density, within a box of side 3/i~^Mpc 



http:/ /www. g- vo.org/Millennium 



centered on each galaxy, and the stellar mass weighted age, 
using the below query. 

SELECT d.galaxylD, d.haloID, d. stellarMass , d.snapnum, s. redshift, 

s . lookBackTirae , dl.galaxylD as galaxyid2 , d .massWeightedAge , 

dl .raassWeightedAge AS massWeightedAgeGal , d, centralMvir ,sdss.r_sdss 

FROM DeLucia2006a d JOIN Snapshots s ON s . snapnum=d. snapnum 
AND d.snapnum=60 JOIN DeLucia2006a dl ON d, snapnuni^dl . snapnum 
AND dl . snapnuiii=60 AND ABS(d . x-dl . x) <1.5 AND abs (d . y-dl . y) <1.B 
AND ABS(d.z-dl .z) <1.5 AND ABS (d. centralMvir - Msval) < BO 
JOIN DeLucia2006a_SDSS2MASS AS sdss on dl . galaxyid=sdss . galaxyid 

We called this query multiple times by looping over in- 
creasing values of "Msval" to obtained all halos with cen- 
tral halo dark matter mass within | d. centralMvir - Msval 
|< 50. This was performed to circumvent timeout errors. We 
chose to obtain data from the simulations with snap numbers 
60, 58, 56, 51, 50, 48 which roughly correspond to the centers 
of the redshift bins used in the analysis which follows. 



3 METHODOLOGY 

In this section we describe the construction of the local 
galaxy densities and the assignment of dark matter halo 
masses to the observational and simulated semi analytic 
data. We note that here the density is only dependent on 
cosmology through the volume element, which changes by 
less than 6% (9%) from 0.2 < z < 0.5 in a fiat fiducial 
ACDM when changing the cosmological parameter Q,rn from 
0.3 to 0.25 (or h from 0.7 to 0.72). 



3.1 Observational local galaxy density 

The magnitude limits of the SDSS imply that we cannot 
directly compare the number density of all observed galaxies 
at different redshifts. We therefore determined the maximum 
volume out to which each galaxy could have been observed, 
and weighted the density by this factor. This is the standard 
Vmax correction ( Schmidt]] 1968 1 and the maximum volume 
was calculated by spline fitting to the faint end of the band 
absolute magnitude Mr and redshift z distributions of the 
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MGS and LRGs separately. The equation of the sphne fit is 
given by 2 = 6.f 2 + 0.91M^ + OMM^ + 7 x W'^M^ for the 
LRGs and z = -262.22 - 31.52M^ - 1.26M^ - O.OieM^^ for 
the MGS. We additionally removed galaxies fainter than the 
spline fit. 

We calculated the local density in a cylindrical cell of 
radius ro = 2.25^~'^Mpc (i.e. diameter 4.5/i~'^Mpc) and 
length lo = 4.5h~^Mpc and additionally weighted each 
galaxy by its Vmax correction. We also examined the effect 
on our results of changes in the cylinder volume, by modify- 
ing both the radius and length to ri = 1.5 &c h = 3 Mpc 
and r2 = 3 &i h = Qh~^ Mpc. We note that the dependence 
of stellar MWA with density for the LRGs using both stel- 
lar population models, are indifferent to the changes in the 
cylinder volumes. For the BC03 model applied to the MGS, 
the stellar MWA at lower densities also remain unchanged, 
but at higher densities there is an increase of ~ 0.5 Gyrs 
for radii increases from ro to r2, and the MWA decreases 
by ~ 1 Gyrs for radii decreases from ro to ri. For the M05 
models applied to the MGS, the stellar MWA also remains 
unchanged at low densities, but at higher densities it de- 
creases by 1.5 Gyrs for radii changes from ro to ri. 

Edge and survey holes were accounted for by creating 
a pixelised mask us ing the Hierarchica l Equal Area isoLati- 
tude Pixelizatior[^ ( |G6rski et al.|2005| hereafter HEALPix) 
routine, applied to a large SDSS DR7 photometric galaxy 
sample with resolution Nside= 128, and only retaining pix- 
els if they were completely surrounded by other pixels. At 
this HEALPix resolution there are many pixels per cylindri- 
cal cell. We divided the galaxy distributions into comoving 
distances slices such that there were three slices per cylinder 
length. 

Each pixel in each redshift slice, which contained at 
least one galaxy became the center of a cylinder (which it- 
self spanned the redshift slice in front and behind of the red- 
shift slice in question), within which we calculated the Vmax 
weighted density for all galaxies in the cylinder (irrespective 
of the number of recovered VESPA age bins), and we fur- 
ther divided this by the total number of galaxies within the 
three redshift slices encompassing the cylinder. We then cal- 
culated the average stellar MWA of only the galaxies which 
resided in the central pixel. 



3.2 Observationally-estimated dark matter halo 
masses 

We matched the SDSS DR7 galaxies in our sample to 



the group members catalogue ( |Yang et al. 20071 of the 



New York University Value-Added Galaxy Catalogue (NYU- 
VAG C|Blanton fc others. |2005|), wh ich was run on the SDSS 
DR4 ( [Adelman-McCarthy et al.|20 06) Main Galaxy Sample. 
The group catalogue prescribes a dark matter halo mass to 
galaxy groups and clusters as a function of the total opti- 
cal luminosity of galaxies above a characteristic luminosity 



(see Yang et al. 2007 for details). This optical-based dark 
matter halo mass was found to agree with the "true" dark 
matter halo mass when compared with mock galaxy redshift 
surveys. For each dark matter halo, we calculated the mean 
of the stellar MWA of galaxy members. 



http:/ /healpix.jpl. nasa.gov 



As we were only able to estimate the dark matter halo 
masses of the low redshift MGS, we have restricted our com- 
parison of the galaxies from the semi analytic models, to the 
low redshift snapshots of the Millennium Simulation. We 
calculated the average stellar MWA of all galaxies in each 
parent dark matter halo. 



3.3 Local galaxy density and halo masses from 
simulations 

Using the SDSS data we determined the maximum abso- 
lute r band magnitude of galaxies within each redshift slice, 
and applied the same absolute magnitudes cuts to the galax- 
ies drawn from the semi analytic models. These cuts were 
— 16.24 for the three lower redshift snapshots, —22.66 for 
the higher redshift snapshots. The magnitude cuts applied 
to the higher redshift simulation outputs result in a lower 
number of recovered galaxies and halos. The output of the 
SQL query returned the number of galaxies within a box 
of side 3/i~'^Mpc centered on each galaxy. We then simply 
calculated the number of galaxies in each box which passed 
the absolute magnitude cuts and divided by the box size. 
We then recorded the stellar MWA of the galaxy at the cen- 
ter of the box. We additionally calculated the average stellar 
MWA of all galaxies (which, again had passed the absolute 
magnitude cuts) in each parent dark matter halo. 



4 RESULTS 

4.1 Local galaxy density 

In Fig. [l] we present the stellar MWA of VESPA galaxies 
using the BC03 and M05 models, as a function of Vmax 
weighted density. We also show the stellar MWA of the semi 
analytic model galaxies, as a function of local galaxy density 
at a range of redshifts chosen to coincide with the observed 
VESPA galaxies. The la dispersion of the stellar MWA is 
indicated using the grey hashed regions, and the error on the 
mean by the black error bars. We show the median relations 
for each redshift by the colored lines. 

Both Fig. [T(ay] and [T(b)] show that the SDSS MGS sam- 
ple probe lower local densities than the higher redshift LRG 



sample, as expected from observations (e.g., Bamford et al. 
[2009 p . 

We find that there is a large spread in stellar MWA as 
a function of l/mai^-weighted density for both stellar popu- 
lation models with the VESPA data, and also find a sim- 
ilar stellar MWA dispersion for the semi analytic models, 
see Fig 1(c) although the density parameter has not been 
Knaa;-corrected because the simulations are not magnitude 
dependent. Weighting the semi analytic model galaxies by 
the Vmax correction could add extra scatter. 

We can try to understand if the large dispersion is ac- 
tually a feature of the data, if it is due to systematic offsets 
in the stellar population models, or the noise in the galaxy 
spectra. Comparing the medians (the colored lines) of the 
stellar MWA of the two different synthetic stellar population 
models run on the almost identical input data sets (the SDSS 
galaxies), we see that there is typically a 1 — 2 Gyr difference 
between them, which is indicative of a systematic offset. Ad- 



ditionally we can examine the results of Tojeiro et al. ( 2010 1, 
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Figure 1. The mass weighted age as a function of Vmax weighted density, for SDSS galaxies. Fig |l(a)| shows the mass weighted age 
reconstructed using the BC03 stellar population models, and Fig |l(b)"] uses the M05 models. Fig |l(c)| shows the mass weighted age of 
galaxies drawn from the Millenium Simulations by applying semi analytic (SA) models to the dark matter halos. The higher redshift 
LRGs and lower redshift MGS are each subdivided into three redshift bins, shown by the color lines and black points. We show the la 
standard deviation by the grey region, and the error on the mean with the black error bars. 



who stacked high signal-to-noise spectra from the SDSS DR7 
LRG sample ( |Abazajian et al.|2009[ ), and compared VESPA 
outputs (e.g., redshift evolution of metallicity weighted age, 
MWA) between the BC03 and M05 models, see their Fig 6. 
They also find a > 20% dispersion between the models, but 
also find that both models need to be modified to include 
additional error terms to fit the observed spectra well. These 



results suggest that the dispersion is at least partly due to 
the uncertainty in the models and systematic differences, in- 
dependent of uncertainties in the data. If we now compare 
the dispersion to that seen in the galaxies drawn from the 
semi analytic model, we conclude that the remaining scatter 
is a combination of being intrinsic to the data and due to 
noise. 
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Figure 2. The stellar MWA and Vmax weighted density relation of the low redshift MGS, subdivided into bins of total stellar mass for 
the BC03 model (left panel) and M05 model (right panel). We use the colored lines to distinguish redshift slices, and like colored lines 
of increasing stellar MWA correspond to bins of increasing total stellar mass. The stellar mass M bins are given in units of Mq, by the 
ranges 10.2 ^ logjg ^4" < 10.4, 10.4 ^ logio M < 10.8, 10.8 ^ logio M < 11.2 and 11.2 ^ logio M < 11-6, which split the sample into 
bins of roughly equal number. The error on the mean is again showed by the black error bars, and the \o dispersion by the grey hashed 
region. 



Examining Fig. 1(a) and 1(b) we see that, as noted 
above, the two different steUar population models agree to 
within ~ 20%, and that in both cases there is a ~ 1 Gyr dif- 
ference between the stellar MWA of LRGs and MGS, which 
is a consequence of an older stellar population, i.e. there 
has been comparatively little or no recent star formation 
in LRGs compared with MGS. We note that a direct com- 
parison of the BC03 and M05 stellar MWA for the 21801 
LRG galaxies with measured stellar MWA and the con- 
straint Nbins > 4 has a large scatter (68% = 3.7 Gyrs) 
and a small systematic offset of 0.6 Gyrs. A possible expla- 
nation of the systematic offset, is that, of the two stellar 
population models used in VESPA, the M05 had a lower 
wavelength resolution than BC03. VESPA re-bins the data 
to the models, so the signal-to-noise is larger in the M05 
models, and thus more age bins are recovered. Also there 
is a decrease in star formation between 0.1 and IGyr using 
BC03, which will mean there will be little stellar mass com- 
piled in the corresponding VESPA age bins, i.e., the stellar 
MWA will be skewed low. 



Continuing to examine Fig. 1(a) and 1(b) we also find 
that the stellar MWA for all galaxies is typically older than 
8 — >■ 9.5 Gyrs in look back time for the MGS, and typically 
older than 9 — >• 10 Gyrs for the LRGs. We note that Fig. 



1(c) is divided by redshift, but does not distinguish between 
LRG and MGS galaxy populations (i.e. they are combined 
here), but do produce the same trends, i.e. that the largest 
fraction of stellar assembly occurred before 8 Gyrs look back 
time. 

For each galaxy sample analysed with the BC03 model, 



we find that the median stellar MWA of the three increas- 
ing redshift slices increases. This suggests that even though 
the majority of star formation occurred at high redshift, 
both the LRG and MGS samples continuously have small 
amounts of stars being assembled, thereby reducing their 
stellar MWA. This effect is less apparent with M05 model, 
with only the two upper redshift slices of the MGS obeying 
this trend. 

This weak trend is however consistent with observations 
of the continual Star Formation Rates (SFR) of both low and 



high mass galaxies across cosmic time (Gonroy & Wechsler 



[2009) INoeske et al]|2007a|b[ although note, that the SFR 
for high redshift massive galaxies is greater than for high 
redshift low mass galaxies). 

All panels of Figs, [l] appear to show no, or little depen- 
dence of stellar MWA on local galaxy density to first order. 
Also note that the observed and simulated trends agree with 
each other within the dispersion. Furthermore, we see no dis- 
cernible trends between the different redshift slices of either 
the MGS or LRGs as a function of density. 

A correlation of very local galaxy densi ty (related to 
projected density of the S**" nearest neighbour. Cooper et al. 



2005 1 and galaxy age estimates following Gallazzi et al 



( 2005 1 have led Cooper et al. ( 2010 1 to find a strong relation 



between galaxy age and local galaxy density at fixed stellar 
mass. Our estimates of density are in much larger volumes, 
and thus these smaller scale correlations are smoothed out. 

In Fig. |2] we have split the stellar MWA and Vmax- 
weighted density relation into bins of total stellar mass 
for both models. We use the colored lines to distinguish 
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(a) Dark matter MAH 




11,5 12.0 13,5 13.0 13,5 14,0 14,5 15,0 11.0 11.5 12,0 12.5 13.0 13,5 14,0 

Group Mass [log^pM^] Halo Mass [log^j, M^] 

(b) BC03 & M05 models, with DM MWA (c) Semi analytic models 

Figure 3. We show the dark matter assembly history, using the MAH code, and the assembly history in age bins reconstructed as could 
be measured with VESPA in Fig |3(a)] Fig |3(b)"] shows the stellar MWA as a function of group halo mass for both the M05 models (red) 
and BC03 models (green), and the la dispersion by the grey hashed region, and the error on the mean by the black lines. We also show 
the dark matter MWA (DM MWA), as calculated by MAH code (black dotted lines) and show the Itr dispersion on the dark matter 
MWA with the blue regions, which are determined by 10* random re-samples using the bin ages available to VESPA. We additionally 
show, in Fig |3(c)| the stellar MWA as a function of halo mass, for galaxies from the semi analytic models in the Millennium Simulation, 
and again show the la dispersion by the grey hashed region, and the error on the mean by the black error bars. 



redshift slices, and like colored lines of increasing stellar 
MWA correspond to bins of increasing total stellar mass. 
The stellar mass M bins are given in units of Mq, by the 
ranges 10.2 ^ logjo M < 10.4, 10.4 logm M < 10.8, 
10.8 sC logioM < 11.2 and 11.2 ^ logn, M < 11.6, which 
split the sample into bins of roughly equal galaxy number. 



The error on the mean is again showed by the black er- 
ror bars, and the la dispersion by the grey hashed region. 
We have only shown the mass splitting of the lower redshift 
MGS to make the figures more readable. There is very little 
stellar MWA splitting between the bins of equal mass for 
each redshift slice of the LRG sample. 
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The stellar MWA is consistent with being independent 
on local galaxy density, independent of total stellar mass 
bin, redshift slice or stellar population model, within the 
la spread of the data. However, we do find that the stellar 
MWA of lower-mass galaxies is consistently lower than more 
massive galaxies, across all redshift slices and independent 
of the stellar population model. This suggests that the lower 
the mass of a galaxy, the more stellar mass is compiled at 
late times, i.e., low mass galaxies are younger than high 
mass galaxies. This result has been observed before, e.g., 



Jimenez et al. (20051; Panter et al. (20071 found that the 



star formation rates of massive galaxies peaked earlier than 
less massive galaxies. 

4.2 Dark Matter Halo Mass 

We now compare both the stellar MWA (for low-redshift 
galaxies matched to the DR4 NYVAC group catalogue) and 
the dark matter MWA, to the dark matter host halo mass. 
To do so we used the publicly available code of |Zhao et al.| 
( 2009 1 which calculates the Mass Accretion Histories (here- 



after MAH) of dark matter halos back through cosmic time, 
by assuming that the progenitor of a halo is always the most 
massive of all progenitors. The fitting formula was calibrated 



on suites of N-body dark matter simulations (see Zhao et al. 
2009[ for details) and was found to be accurate to within 



< 5% and universal. Results may be defined in terms of 
M200, which enabled direct comparison with the halo masses 
of the NY VAC group catalogue. 

In Fig. 3(a) we show, as a function of halo mass, the 



percentage of accreted dark matter onto the dark matter 
halos per time interval, relative to the halo mass at 2; = 0. 
The black lines show the time steps as outputted by the 
MAH code, and the red lines show different realizations of 
the larger time steps corresponding to a random selection 
of > 4 VESPA age bins. The VESPA age bins were cho- 
sen by imposing the constraint that each random sample of 
bin selections be consecutive and complete across the full 
VESPA time range. This age bin selection was performed 
to most closely match the conditions of the VESPA data; 
recall that each included galaxy spectrum has a high signal- 
to-noise, such that VESPA could determine stellar histories 
in > 4 different age bins. The different line styles correspond 
to different mass halos at 2; = 0. 



In Fig. 3(b) we show the average stellar MWA of halo 



member galaxies as a function the dark matter group halo 
mass, and distinguish between the M05 models (red) and 
BC03 models (green). We show the la dispersion by the grey 
hashed region, and the error on the mean by the black error 
bars. We again note that the spread of data is large, and 
that the models agree to within ~ 1.5 Gyrs. For both mod- 
els we see a slight increase in the stellar MWA for increasing 
dark matter halo mass, which suggests that more massive 
groups and clusters formed the majority of their stars earlier 
in the Universe's history than less massive groups, or sim- 
ilarly, that stellar mass was still being assembled in lower 
mass groups more recently. We additionally show the dark 
matter MWA by the black dotted line. The blue contours 
correspond to 10'' dark matter MWA measurements drawn 
from random samplings of the 30 VESPA age bins, with the 
constraints (as introduced above), that we randomly choose 
more than 4 VESPA bins, and that the bins form a full and 



continuous coverage of VESPA ages. We see that the dark 
matter MWA becomes lower as the group mass increases, 
and that the la simulated dispersion is constant as a func- 
tion of group mass. 

The stellar and dark matter MWA of the smaller sys- 
tems have similar values, but this may not be a physical 
correlation because of the anti-correlation of the stellar and 
dark matter MWA at higher masses. This leads us to infer, 
that while the stellar MWA is correlated with the overall 
mass of the dark matter halo, is it not correlated with the 
accretion of dark matter onto the halo. 



Fig 3(c) shows the stellar MWA as a function of dark 
matter halo mass, as measured from galaxies using semi ana- 
lytic models in the Millennium Simulation. We note that the 
masses probed by the simulations is lower than those from 
the NYVAC, but we do see a clear increase in the stellar 
MWA with increasing dark matter halo mass, in agreement 
with the VESPA data, and again opposite to the trends seen 
in the dark matter MWA. 



5 CONCLUSIONS 



We used the results of VESPA (iTojeiro et al."2007! '2009 
which analyzed 10^ SDSS DR7 (Abazajian ct al. 200j 
galaxy spectra, to calculate the stellar Mass Weighted Age 
(MWA). To remove peculiarities in the data, we chose to 
only use high signal-to-noise spectra which allowed the mea- 
surement of stellar histories in greater than 4 VESPA age 
bins. We then calculated the local galaxy density of the all 
the SDSS galaxies in cylinders of radius 2.25/i~'^Mpc and 
length 4.5h~^ Mpc, and applied the standard Vmax correc- 
tion to account for the magnitude limits of the SDSS. We 
also cross-matched the VESPA galaxies with the New York 



Value Added Group Catalog (Yang et al. 20071, which in- 



cludes estimates of the parent dark matter halo mass. 



Using the public code of Zhao et al. (2009), we deter- 
mined the dark matter MWA as a function of halo mass. We 
additionally obtained the stellar MWA, local galaxy den- 
sity, and dark matter halo masses of galaxies from within 
the Millennium Simulation ( Springel et al.|2005[), populated 
according to the semi analytic models of |De Lucia et al.| 
(2006). 

We found a large scatter in the recovered stellar MWA, 
which we determined to be partially due (~ 20%) to the dif- 
ferences between the BC03 ( [Bruzual fc Charlot||2003| and 
M05 (Maraston et al. 20091 stellar populations models, by 



comparing the medians of our distributions, and also by ex- 
amining the results of [Tojeiro et al.| ( |2010[ ), who compared 
the differences between models when applied to stacked LRG 
spectra with high signal-to-noise. A large dispersion is also 
observed in the stellar MWA of the semi-analytic models, 
suggesting that the dispersion is also inherent to the data. 
We found that the recovered stellar MWA using the different 
stellar population models agree to within 1.5Gyrs and that 
the stellar MWA of most (60% to 90% depending on stel- 
lar population model and galaxy sample) galaxies was older 
than 8 Gyrs, independent of redshift, local galaxy density, 
dark matter halo mass or galaxy type. This is expected from 
the star formation history of the Universe, which peaks at 
around z = 1 — 2, corresponding to 8 — 10 Gyrs lookback 
time (e.g., see Heavens et al.|2004 i. 
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We now return to the questions posed in the introduc- 
tion; 

(i) Does the stellar MWA depend on the local density? 

(ii) Are similar trends seen in semi-analytic models? 

To first order, we found that stellar MWA does not appear 
to be related to local galaxy density in either the observed or 
simulated data. We found similar dispersions in the observed 
and simulated data and apparent flatness across local galaxy 
density (although see Cooper et al.|2010 1 . 



(iii) Is the stellar MWA correlated with the mass of the 
Dark Matter Halo that the galaxy inhabits? 

We did find a correlation of older stellar MWA with increas- 
ing dark matter halo mass in the observed galaxy sample, 
independent of stellar population model, which was also seen 
in the semi-analytic models, albeit for a smaller range in 
dark matter halo mass. 

(iv) Is the stellar MWA correlated with the dark matter 
MWA? 

We found that the dark matter MWA became anti correlated 
with the stellar MWA as the mass of dark matter halos 
increases. This is an observation of a "downsizing" effect 
( Cowie et al.|[T996 l which describes the idea that the dark 
matter halo mass at which star formation is highest, shifts 
to lower masses at later times (see e.g., Fig 



9 of Conroy 

Wechsler|2009| [Jimenez et al.|2005| ) . Another meaning of 



"downsizing" is that more massive galaxies compiled more 
of their stars at higher redshift (and in a shorter time scale) 
than less massive galaxies as seen in Fig. [2] and e.g., [Thomas] 
et al. (20051; Jimenez et al. (20071. These results have been 



seen previously by measuring the Specific Star Formation 
Rates (the amount of star formation per solar mass) of sets 
of similar mass galaxies over a range of redshifts. 

Here, the downsizing effect is observed by the decrease 
in the stellar MWA as a function of decreasing stellar and 
dark matter halo mass, using VESPA stellar history re- 
constructions of the galaxy spectra. These conclusions are 
in agreement with other measures of "downsizing" using 
MOPED to reconstruct the Star Formation Rate (SFR) his- 
tories of SDSS DR3 galaxies ( [Panter et al]|2007| found the 
SFR of massive galaxies peaked earlier than less massive 
galaxies). 

We find VESPA is well suited for reconstructing stel- 
lar histories using the included synthetic stellar population 
models, and agrees well with semi-analytical models drawn 
from larger simulations. 



ACKNOWLEDGMENTS 

The authors would like to thank Rita Tojeiro for VESPA 
database support and useful suggestions, and Aday Robaina 
for insightful discussions and feedback which have improved 
the manuscript. BH would like to thank the Theory Divi- 
sion, Department of Physical Sciences, University of Helsinki 
for office space, and acknowledges funding from grant num- 
ber FP7-PE0PLE- 2007- 4-3-IRG n 20218, LV is supported 
by FP7-PEOPLE-2007-4-3-IRG n. 202182, FP7-IDEAS- 
Phys.LSS 240117; LV and RJ are supported by MICINN 



grant AYA2008-03531 Funding for the SDSS and SDSS- 
II has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Science Foun- 
dation, the U.S. Department of Energy, the National Aero- 
nautics and Space Administration, the Japanese Monbuk- 
agakusho, the Max Planck Society, and the Higher Educa- 
tion Funding Council for England. The SDSS Web Site is 
http://www.sdss. org/[ 



REFERENCES 

Abazajian K., et al., 2005, AJ, 129, 1755 
Abazajian K. N., et al., 2009, ApJS, 182, 543 
Adelman-McCarthy J. K., et al., 2006, ApJS, 162, 38 
Bamford S. P., et al., 2009, MNRAS, 393, 1324 
Blanton M. R., others. 2005, AJ, 129, 2562 
Boylan-Kolchin M., Springel V., White S. D. M., Jenkins 

A., Lemson G., 2009, MNRAS, 398, 1150 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Conroy C., Wechsler R. H., 2009, ApJ, 696, 620 
Cooper M. C, Gallazzi A., Newman J. A., Yan R., 2010, 

MNRAS, 402, 1942 
Cooper M. C., Newman J. A., Madgwick D. S., Gerke B. F., 

Yan R., Davis M., 2005, ApJ, 634, 833 
Cowie L. L., Songaila A., Hu E. M., Cohen J. G., 1996, AJ, 

112, 839 

Croton D. J., et al., 2006, MNRAS, 365, 11 
De Lucia G., Blaizot J., 2007, MNRAS, 375, 2 
De Lucia G., Kauffmann G., White S. D. M., 2004, MN- 
RAS, 349, 1101 
De Lucia G., Springel V., White S. D. M., Croton D., Kauff- 
mann G., 2006, MNRAS, 366, 499 
Eisenstein D. J., et al., 2001, AJ, 122, 2267 
Ferreras I., Pasquali A., Rogers B., 2010, ArXiv e-prints 
Gallazzi A., Chariot S., Brinchmann J., White S. D. M., 
Tremonti C. A., 2005, MNRAS, 362, 41 
Gorski K. M., et al., 2005, ApJ, 622, 759 
Graff P., Hobson M., Lasenby A., 2010, ArXiv e-prints 
Heavens A., Panter B., Jimenez R., Dunlop J., 2004, Na- 
ture, 428, 625 

Heavens A. F., Jimenez R., Lahav O., 2000, MNRAS, 317, 
965 

Jimenez R., Bernardi M., Haiman Z., Panter B., Heavens 
A. F., 2007, ApJ, 669, 947 

Jimenez R., Panter B., Heavens A. F., Verde L., 2005, MN- 
RAS, 356, 495 

Maraston C, 2005, MNRAS, 362, 799 

Maraston C, Stromback G., Thomas D., Wake D. A., 
Nichol R. C., 2009, MNRAS, 394, L107 

Mateus A., Jimenez R., Gaztanaga E., 2008, ApJL, 684, 
L61 

Noeske K. G., et al., 2007a, ApJL, 660, L47 
Noeske K. G., et al., 2007b, ApJL, 660, L43 
Panter B., Jimenez R., Heavens A. F., Chariot S., 2007, 

MNRAS, 378, 1550 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 

525 

Schmidt M., 1968, ApJ, 151, 393 
Springel V., 2005, MNRAS, 364, 1105 
Springel V., et al., 2005, Nature, 435, 629 



10 Hoyle et al. 



Springel V., Yoshida N., White S. D. M., 2001, New A, 6, 
79 

Strauss M. A., et al., 2002, AJ, 124, 1810 

Thomas D., Maraston C, Bender R., Mendes de Oliveira 

C, 2005, ApJ, 621, 673 
Tojeiro R., Heavens A. F., Jimenez R., Panter B., 2007, 

MNRAS, 381, 1252 
Tojeiro R., Percival W. J., Heavens A. F., Jimenez R., 2010, 

ArXiv e-prints 



Tojeiro R., Wilkins S., Heavens A. F., Panter B., Jimenez 

R., 2009, ApJS, 185, 1 
Yang X., Mo H. J., van den Bosch F. C, Pasquah A., Li 

C, Barden M., 2007, ApJ, 671, 153 
York D. G., ct al., 2000, AJ, 120, 1579 
Zhao D. H., Jing Y. P., Mo H. J., Borner G., 2009, ApJ, 

707, 354 



